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CN . Abstract 

Q , A surprisingly large asymmetry between the up and down sea quark distributions in the nucleon 

C/^ ■ has been observed in recent deep inelastic scattering and Drell-Yan experiments. This result 

strongly suggests that the mesonic degrees of freedom play an important role in the description of 
the parton distributions of the hadronic sea. In this article, we review the current status of our 
knowledge of the flavor structure of the nucleon sea. The implications of various theoretical models 
J> ■ as well as possible future measurements are also discussed. 
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§s 1 Introduction 



The first direct evidence for point-like constituents in the nucleons came from the observation of scaling 
phenomenon in Deep-Inelastic Scattering (DIS) experiments []T], |^ at SLAG. These point-like charged 
1) I constituents, called partons, were found to be spin-1/2 fermions. These partons were initially identified 
as the quarks in the Constituent Quark Models (CQM). However, it was soon realized that valence 
^ [ quarks alone could not account for the large enhancement of the cross sections at small Bjorken-x, 
^ I the fraction of nucleon momentum carried by the partons. This swarm of low-momentum partons, 
dubbed wee-partons by Feynman 0, was interpreted as the quark and antiquark sea of the nucleon. 
^ . DIS experiments therefore provided the first evidence for the existence of antiquarks in the nucleon. 

The observation of partons in DIS experiments paved the road to the formulation of Quantum 
Chromodynamics (QCD) as the theory for strong interactions. Nevertheless, the exact form of the 
parton distribution functions (PDF) can not be deduced from perturbative QCD. Like many other 
static properties of hadrons, the parton distribution functions belong to the domain of non-perturbative 
QCD. In spite of great progress made in Lattice Gauge Theory (LGT) in treating the bound-state 
properties of hadrons, it remains a challenge to predict the parton distributions using LGT. 

Until parton distributions can be readily calculated from first principles, they are best determined 
from experiments. Electroweak processes such as DIS and lepton-pair production provide the cleanest 
means to extract information on the parton distributions. There are at least two reasons why it is 
important to measure the parton distribution functions. First, the description of hard processes in 
high energy interactions requires parton distribution functions as an essential input. Second, many 
aspects of the parton distributions, such as sum rules, scaling-violation, asymptotic behaviors at large 
and small x, and flavor and spin structures, can be compared with the predictions of perturbative as 
well as non-perturbative QCD. 

In this article, we review the status of our current knowledge of the flavor dependence of the sea 
quark distributions in hadrons. The recent observation of a striking flavor asymmetry of the nucleon 
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sea has profound implications on the importance of meson degrees of freedom for the description of 
parton substructures in the nucleon. In Section 2, we review the early studies of the nucleon sea in DIS 
and lepton-pair production. The crucial recent experiments establishing the up/down flavor asymmetry 
of the nucleon sea are discussed in Section 3. The various theoretical models for explaining the flavor 
asymmetry of the nucleon sea are described in Section 4. The implications of these models on other 
aspects of the parton structure functions are discussed in Section 5. Finally, we present future prospects 
in Section 6, followed by the conclusion in Section 7. 



2 Early Studies of the Nucleon Sea 

2.1 Deep Inelastic Scattering 

Although scaling behavior in inelastic electron scattering was predicted by Bjorken based on the 
framework of current algebra, its confirmation by the SLAC experiments still came as a major surprise. 
A simple and intuitive picture for explaining the scaling behavior is the parton model advanced by 
Feynman P, P]. In this model, the electron- nucleon deep- inelastic scattering is described as an incoherent 
sum of elastic electron-parton scattering. However, the nature of the partons within the nucleon was not 
specified by Feynman. Some authors 0, H, |T0| speculated that the partons were the 'bare nucleon' plus 



the pion cloud, while others [|ri|, |T^, believed they were the quarks introduced by Gell-Mann 
and Zweig [0. The latter scenario was strongly supported by the measurement of R, ratio of the 
longitudinally over transversely polarized photon cross sections, showing the spin-1/2 character of the 
partons. 

Evidence for quark-antiquark sea in the nucleon came from the observation that the structure 
function F2{x) approaches a constant value as x — |T6[. If the proton is made up of only three 



quarks, or any finite number of quarks, F2{x) is expected to vanish as x ^ 0. Bjorken and Paschos 
therefore assumed that the nucleon consists of three quarks in a background of an infinite number 
of quark-antiquark pairs. Kuti and Weisskopf |jl2| further included gluons among the constituents of 



nucleons in order to account for the missing momentum not carried by the quarks and antiquarks alone. 

The importance of the quark-antiquark pairs in the nucleon is in sharp contrast to the situation for 
the atomic system, where particle-antiparticle pairs play a relatively minor role (such as the polarization 
of the vacuum). In strong interactions, quark-antiquark pairs are readily produced as a result of the 
relatively large magnitude of the coupling constant as, and they form an integral part of the nucleon's 
structure. 

Neutrino-induced DIS experiments allowed the separation of sea quarks from the valence quarks. 
Recall that 

Fr(x) = 2a;^[g,(x)+g,(x)], 



where i denotes the flavor of the quarks. Note that the valence quark distribution is defined as the 
difference of the quark and antiquark distributions, = qi{x) — qi{x). Equation |I] shows that the 

valence quark distribution is simply F^^{x)/2, while the sea quark distribution is given by F2^{x)/2x — 



F!^^{x)/2. The F2^{x) and F^^{x) data from the CDHS experiment [17| clearly showed that the valence 
quark distributions dominate at x > 0.2, while the sea quarks are at small x. 

The earliest parton models assumed that the proton sea was flavor symmetric, even though the 
valence quark distributions are clearly flavor asymmetric. Inherent in this assumption is that the 
content of the sea is independent of the valence quark's composition. Therefore, the proton and neutron 
were expected to have identical sea-quark distributions. The assumption of flavor symmetry was not 
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based on any known physics, and it remained to be tested by experiments. Neutrino-induced charm 
production experiments [|18|, |l9l, which are sensitive to the s ^ c process, provided strong evidences 
that the strange-quark content of the nucleon is only about half of the up or down sea quarks. This 
flavor asymmetry was attributed to the much heavier mass for strange quark compared to the up and 
down quarks. The mass for the up and down quarks being very similar suggests that the nucleon sea 
should be nearly up-down symmetric. A direct method to check this assumption is to compare the sea 
in the neutron to that in the proton by measuring the Gottfried integral in DIS, as discussed next. 



2.2 Gottfried Sum Rule 

In 1967, Gottfried studied electron-proton scattering with the assumption that the proton consists 
of three constituent quarks He showed that the total electron-proton cross section (elastic plus 



inelastic) is identical to the Rutherford scattering from a point charge. Gottfried derived a sum rule 

/f = Q')/x dx = Y.m' = 1, (2) 

Jo . 

where is the charge of the ith quark in the proton. Gottfried expressed great skepticism that this 
sum rule would be confirmed by the forthcoming SLAG experiment by stating "/ think Prof. Bjorken 
and I constructed the sum rules in the hope of destroying the quark modeT |^ . Indeed, Eq. ^ was 
not confirmed by the experiments, not because of the failure of the quark model, but because of the 
presence of quark- ant iquark sea. In fact, the total number of the sea partons being infinite makes /f 
diverge. A closely related sum rule, now called the Gottfried Sum Rule (GSR), avoids this problem by 
considering the difference of the proton and neutron cross sections, namely, 

II - ^ = fmx^Q') - F-{x,Q^)\ix dx = j:[m' - (gn^i = vs. (s) 

Jo Y 

In deriving Eq. |^, it was assumed that the sea quarks in the proton and neutron are identical. Regarding 
the Gottfried Sum Rule, Kuti and Weisskopf stated in their paper |jl2[ ''This very simple and definite 
sum rule is based upon our most radical simplification: the assumption that the core carries vacuum 
quantum numbers so that the isotopic spin of the nucleon is completely carried by the valence quarks" . 

Soon after the discovery of scaling in electron-proton DIS, electron-deuterium scattering experiments 
were carried out to extract the electron-neutron cross sections. The comparison of e—p with e — n data 



was very important for distinguishing early competing theoretical models [16|. These data also allowed 



a first evaluation ||2^ of the Gottfried integral in 1970. The first result for the Gottfried integral was 
0.19, considerably less than 1/3. As the data only covered x > 0.08, it was assumed that Ff — -F^ 
follows Regge behavior (proportional to x^/^) in the unmeasured small-x region. Due to the 1/x factor 
in the integrand, the small-x region could have potentially large contributions to the Gottfried integral. 
Moreover, it was not clear if — would indeed follow the Regge behavior at small x, and if so, at 
what value of x would it set in. By 1973, new data were available down to x = 0.05 and the Gottfried 
integral was evaluated to be 0.28 |2^, considerably larger than the first result. It should be pointed out 



that these data were taken at relatively low values of Q^. Furthermore, varied as a function of x. 

Although the large systematic errors associated with the unmeasured small-x region prevented a 
sensitive test of the GSR, Field and Feynman nevertheless interpreted the early SLAG data as a 



strong indication that GSR is violated and that the u and d distributions in the proton are different. 
The relationship between the Gottfried integral and the d/u asymmetry is clearly seen in the parton 
model, namely. 



f\F^{x, Q^) - F^{x, Q'')]/x dx = l + l f\u{x, Q^) - J(x, Q^)]dx. (4) 
o 3 JO 
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Figure 1: EMC El and SLAC pg measurements of - F^. 



Equation ^ clearly shows that the early SLAC data on GSR implied d > u, a.t least for certain region 
of X. Field and Feynman further suggested that Pauli blocking from the valence quarks would inhibit 
the uu sea more than the dd sea, hence creating an asymmetric nucleon sea. 

The SLAC DIS experiments were followed by several muon-induced DIS experiments at Fermilab 
and at CERN. Using high energy muon beams, these experiments reached much larger values of and 
they clearly observed 26, 27] the scaling- violation phenomenon in DIS. The Gottfried integral was 
also evaluated in muon DIS experiments ^8|, Figure ^ compares the data from the European Muon 
Collaboration (EMC) ||28| with earlier electron data from SLAC The coverages in x are similar 



in these two experiments, even though the values covered by the EMC are much larger. Figure |I| 
shows that F2 — F2 from EMC tend to shift towards smaller x relative to the SLAC data, in qualitative 
agreement with the QCD Q^-evolution. The Gottfried integral determined from the EMC experiment 
is 0.235 + 0.110 — 0.099, consistent with the result from SLAC, but still lower than 1/3. 

Despite the fact that all measurements of the Gottfried integral consistently showed a value lower 
than 1/3, the large systematic errors prevented a definitive conclusion. As a result, all parametriza- 
tions pTI, B2|, [33, p3, BH] of the parton distributions based on global fits to existing data before 1990 



assumed a symmetric u, d sea. As discussed next, more compelling evidences for an asymmetric light- 
quark sea were provided by results from new experiments. 



2.3 Drell-Yan Process 

The first high-mass dilepton production experiment |3^ was carried out at the AGS in 1969, soon after 
scaling was discovered at SLAC. Drell and Yan interpreted the data within the parton model, in 
which a quark-antiquark pair annihilate into a virtual photon subsequently decaying into a lepton pair. 
This simple model was capable of explaining several pertinent features of the data, including the overall 
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magnitude of the cross sections, the scahng behavior of the cross sections, and the polarization of the 
virtual photons. The high-mass continuum lepton-pair production is therefore called the Drell-Yan 
(DY) process. 

Since the underlying mechanism for the DY process involves the annihilation of a quark with an 
antiquark, it is not surprising that this process can be used to probe the antiquark contents of the beam 
or target hadrons. In the parton model, the DY cross section is given by 

Z!eaba(a;i)ga(a;2) +ga(a;i)ga(a;2)]. (5) 



dW^dxF 9M^s (xi + X2) 

Here qa{x) and qa{x) are the quark and antiquark parton distributions of the two colliding hadrons 
evaluated at the momentum fraction x. The sum is over quark flavors. In addition, one has the 
kinematic relations, 

r = X1X2 = M'^/s, 

Xf = Xi ~ X2, (6) 

where M is the invariant mass of the lepton pair and s is the square of the center-of-mass energy. The 
cross section is proportional to a^, indicating its electromagnetic character. 

Equation |^ shows that the antiquark distribution enters as a multiplicative term in the DY cross 
section rather than an additive term in the DIS cross section. Hence, the antiquark distributions can be 
sensitively determined in the DY experiments. The dimuon data from the FNAL E288, in which the T 



resonances were discovered [^, were analysed to extract the sea quark distribution in the nucleon. 
By assuming a flavor-symmetric nucleon sea, namely, u{x) = d{x) = s{x) = Sea{x), the dimuon mass 
distribution obtained in 400 GeV proton-nucleus interaction was described by xSea{x) = 0.6(1 — a;)^°. 



In a later analysis taking into account additional data at 200 and 300 GeV, the E288 collaboration 



found that a much better fit could be obtained with an asymmetric sea, namely, 

u{x) = (1 - s{x) = {u{x) + d{x))/4. (7) 



The need for an asymmetric u and d was also revealed in the E288 da/dy DY data at ?/ = 0, where 



y is the center-of-mass rapidity. For p + A collision, the slope of da/dy at y = is expected to be 
positive due to the excess of u over d valence quarks in the proton. The E288 data showed that the 
slopes are indeed positive, but larger than expected from a flavor symmetric sea. A surplus of d over u 
in the proton sea would lead to more positive slope in agreement with the data pUf . 

The FNAL E439 collaboration |^T| studied high mass dimuons produced in p + W interaction at 
400 GeV. Their spectrometer covered a considerably larger range in xp than E288. They again found 
that an asymmetric sea, u{x) = (1 — x^'^dlx), could well describe their data. 

With all the tantalizing evidence for an asymmetric sea from DIS and DY experiments, it is curious 



that all global analyses pT| , p^ , |33| , |3^ , |35| of parton distributions in the 1980 's still assumed a symmetric 
light-quark sea. This probably reflected the reluctance to adopt an unconventional description of the 
nucleon sea without compelling and unambiguous experimental evidence. As discussed in the next 
Section, such evidence became available in the 1990 's. 

3 Recent Experimental Developments 

3.1 NMC Measurements of the Gottfried Integral 



After the discovery of the so-called 'EMC effect' ||4^ which showed that the parton distributions in 
heavy nuclei are different from that in the deuteron, the EMC detector system was modified by the 
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New Muon Collaboration (NMC) at CERN to study in detail the EMC effect. Special emphases were 
placed on the capability to reach the low-x region where the 'shadowing effect' |^3[ is important, and 
to measure cross section ratios accurately By implementing a 'small angle' trigger which extended 
the scattering angle coverage down to 5 mrad, the lowest value of x reached by NMC was ~ 0.001. The 
NMC also placed two targets in the muon beam, allowing DIS data from two different targets to be 
recorded simultaneously, thus greatly reducing the beam flux normalization uncertainty. To account for 
the different geometric acceptances for events originating from the two targets, half of the data were 
taken using a target configuration where the locations of the two targets were interchanged. 

The upgraded NMC detectors allowed a definitive study |^ of the shadowing effect at low x. 
Moreover, they enabled a much more accurate determination of the Gottfried integral. Figure shows 

in which the smallest x reached (0.004) was significantly 



the Fi - reported by the NMC in 1991 [|6 



lower than in previous experiments. Taking advantage of their accurate measurements of the F^ / F2 
ratios, NMC used the following expression to evaluate — F2-, namely. 



F! - F^ 



F^ 



l-F2VF^)/{l + F2VF!). 

The ratio F^/F^' = F^/Ff — 1 was determined from NMC's F^/Ff measurement, while F2 {Fr 
F2) was taken from a fit to previous DIS experiments at SLAC |^ , Fermilab 
The value of the Gottfried integral for the measured region at = 4 GeV 



(8) 



27 



and CERN ||, H. 
is ^g(0.004 - 0.8) = 



0.227±0.007(stat)±0.014(s?/st). The contribution to from x > 0.8 was estimated to be 0.002±0.001. 
Assuming that F^ - F^ at x < 0.004 behaves as ax\ NMC estimated SciO - 0.004) = 0.011 ± 0.003. 
Summing the contributions from all regions of x, NMC obtained Sq = 0.240 ± 0.016, which was 
significantly below 1/3. This represented the best evidence thus far that the GSR was indeed violated. 
In 1994, NMC reevaluated |^| the Gottfried integral using a new F2 parametrization from Ref. |51 



and newly determined values of F^/Ff. The new F^ parametrization include data from SLAC 
BCDMS I^H, as well as NMC's own measurement ^1|]. The new NMC results for Ff — F^ are shown in 
Fig. ^. Note that the 1994 values are slightly larger (smaller) than the 1991 values at small (large) x. 
The new evaluation gave 5g(0.004 - 0.8) = 0.221 ±0.008(stot) ±0.019(s|/st), and the Gottfried integral 
became 

Sg = 0.235 ± 0.026. (9) 

This value is consistent with the earlier number. The new systematic error is larger than the old one, 
reducing somewhat the overall significance of the NMC measurement. Nevertheless, the violation of the 
GSR is established at a ia level. 



More recently, NMC published their final analysis of F^/Ff [|5^ and F2 ^ 



This analysis 

included the 90 and 280 GeV data taken in 1986 and 1987, as well as the 1989 data at 120, 200 and 280 
GeV. The 1989 data were not used in the earlier evaluations Q of the Gottfried integral. Based 
on these new data, NMC reported ^g(0.004 - 0.8) = 0.2281 ± 0.0065(stat) at = 4 QeV^ H. This 
agrees within statistical errors with previous NMC results [ffS, |5( 



QCD corrections to various parton-model sum rules have been reviewed recently |^5|. The as and 
corrections to the Gottfried integral have been calculated [^, |^ and found to be very small (roughly 
0.4 % each at = 4 GeV^). Therefore, QCD corrections can not account for the large violation of 
GSR found by the NMC. Although perturbative QCD predicts a weak dependence for the Gottfried 
integral, it has been suggested |5^, ^ that due to the non-perturbative origin of the d, u asymmetry 
the dependence of the Gottfried integral will be anomalous between 1 and 5 GeV^. This interesting 
suggestion remains to be tested by DIS experiments. 



3.2 E772 Drell-Yan Experiment 

The main goal of the Fermilab experiment E772 was to examine the origin of the EMC effect. Among 
the many theoretical models which explain the EMC effect |6^, the pion-excess model |61| predicted a 
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Figure 2: NMC measurements of — FJ* (solid data points) and the Gottfried integral 
(open data points). The triangular data points correspond to results published in 1991 [|^, 
while the circular data points represent a more recent analysis in 1994 The extrapolated 
value of Gottfried integral (So) and the expected GSR value are also indicated. 



large enhancement of antiquark content due to the presence of additional meson cloud in heavy nuclei. 
This prediction could be readily tested by measuring the nuclear dependence of proton-induced DY 
cross sections. Using an 800 GeV proton beam, the DY cross section ratios of a number of nuclear 
targets (C, Ca, Fe, W) were measured relative to deuterium with high statistics over the region 



0.04 < X2 < 0.35. The enhancement of antiquark contents predicted by the pion-excess model was not 



observed, and the E772 results were in good agreement with the prediction of the rescaling model [33 



Information on the d/u asymmetry has also been extracted from the E772 data [Q. At xp > 0.1, 
the dominant contribution to the proton-induced DY cross section comes from the annihilation of u 
quark in the projectile with the u quark in the target nucleus. It follows that the DY cross section (per 
nucleon) ratio of a non-isoscalar target (like ^H, Fe, W) over an isoscalar target (like "^H,^^ C,'^" Ca) is 
given as 

RAix2) = (rAix2)/aisix2) ^1 + [{N- Z)/A][{1 - uix^) / d{x2)) / (1 + n(a;2)Mx2))], (10) 

where X2 is the Bjorken-x of the target partons, IS stands for isoscalar, and N, Z and A refer to the 
non-isoscalar target. The [N — Z)/A factor in Eq. [1^ shows that the largest sensitivity to the d/u could 
be obtained with a measurement of a{p + p)/o'{p + d). Nevertheless, for W target {N — Z)/A = 0.183 
and the E772 aw/c^is data could be used to study the d,u asymmetry. 

Figure |^ shows the E772 DY cross section ratios from the neutron-rich W target over the isoscalar 
targets, and ^^C. Corrections have been applied to the two data points at X2 < 0.1 to account for 
the nuclear shadowing effects in C and W . The E772 data provided some useful constraints on the 



d/u asymmetry. In particular, some early parametrizations ^ of large d/u asymmetry were ruled 
out. Despite the relatively large error bars, Figure ^ shows R > 1.0 at X2 > 0.15, which is consistent 
with d > u in this region. 
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Figure 3: The ratios of W over isoscalar targets DY cross sections from E772 [^]. The dashed 
curve is a calculation using the d/u asymmetric parton distributions suggested in Ref. [ |65[] . 
The solid curve corresponds to a calculation using the MRST distribution functions. 



E772 collaboration also presented the DY differential cross sections foip + d at (M^^) = 8.15 GeV. 
As shown in Fig. ^, the DY cross sections near xp = are sensitive to d/u and the data disfavor a large 
d/u asymmetry. Figure ^ also shows that a recent parton distribution function, MRST which has 
modest d/u asymmetry, is capable of describing the p + d differential cross sections well (see Section 
3.6). 

While the E772 data provide some useful constraints on the values of d/u, it is clear that a measure- 
ment of crDY{p + d)/aDY{p + p) is highly desirable. The {N — Z)/A factor is now equal to —1, indicating 
a large improvement in the sensitivity to d/u. Moreover, the uncertainty arising from nuclear effects 
would be much reduced. It has also been pointed out that d/u asymmetry in the nucleon could be 
significantly modified in heavy nuclei through recombination of partons from different nucleons. Clearly, 
an ideal approach is to first determine d/u in the nucleon before extracting the d/u information in heavy 
nuclei. 



3.3 NA51 Drell-Yan Experiment 

Following the suggestion of Ellis and Stirling [Q, the NA51 collaboration at CERN carried out the first 



dedicated dimuon production experiment to study the flavor structure of the nucleon sea |£9[. Using a 



450 GeV proton beam, roughly 2800 and 3000 dimuon events with M^^ > 4.3 GeV have been recorded, 
respectively, for p + p and p + d interaction. The spectrometer setting covers the kinematic region near 
?/ = 0. At y = 0, the asymmetry parameter, Ady, is given as 



Ady = 



a 



pp 



(4Ay - 1) (A, - 1) + (Ay - 1) (4A, - 1) 



aPP + aP^ (4Ay + 1) (A, + 1) + (Ay + 1) (4A, + 1) 



8 



10 



U I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I u 



800 GeV p+ H 



> 

CD 
O 

T3 




-0.2 -0.1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Figure 4: The p+^H DY differential cross sections from E772 [Q. The dotted (dashed) curve 
is a calculation using the parton distributions from Ref. with (without) d/u asymmetry. 
The sohd curve uses the MRST distribution functions. These are leading-order calculations 
normalized to the data. 



where Ay = uv/dy and Xg = u/d. In deriving Eq. 0, the negligible sea-sea annihilation was ignored 
and the validity of charge symmetry was assumed. At x = 0.18, Ay ~ 2 and according to Eq. |TT| 
^DY = 0.09 for a symmetric sea, A^ = 1. For an asymmetric d> u sea, A£,y would be less than 0.09. 

From the DY cross section ratio, (jPp/o-p'^, NA51 obtained A^y = -0.09 ± Q.Q2{stat) ± 0.025{syst). 
This then led to a determination oiu/d = 0.51 ±0.04(stat) ± (syst) at x = 0.18 and (M^^) = 5.22 GeV. 
This important result established the asymmetry of the quark sea at a single value of x. What remained 
to be done was to map out the x-dependence of this asymmetry. This was subsequently carried out by 
the Fermilab E866/NuSea collaboration, as discused next. 



3.4 E866 Drell-Yan Experiment 

At Fermilab, a DY experiment (E866/NuSea) aimed at a higher statistical accuracy with a much wider 
kinematic coverage than the NA51 experiment was recently completed [^. This experiment measured 
the DY muon pairs from 800 GeV proton interacting with liquid deuterium and hydrogen targets. A 
proton beam with up to 2 x 10^^ protons per 20 s spill bombarded one of three identical 50.8-cm long 
cylindrical target flasks containing either liquid hydrogen, liquid deuterium or vacuum. The targets 
alternated every few beam spills in order to minimize time-dependent systematic effects. The dimuons 
accepted by a 3-dipole magnet spectrometer were detected by four tracking stations. An integrated flux 
of 1.3 X 10^^ protons was delivered for this measurement. 

Over 330,000 DY events were recorded in E866, using three different spectrometer settings which 
covered the regions of low, intermediate and high mass muon pairs. The data presented here are from 
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Figure 5: The ratio a^'^ /2a^^ of Drell-Yan cross sections versus X2 for E866. The curves are 
next-to-leading order calculations, weighted by acceptance, of the Drell-Yan cross section 
ratio using various parton distributions. The CTEQ4M and MRS(R2) parton distributions 
existed prior to the E866 data, while the other distributions came after the E866 results 
were obtained. In the lower CTEQ5M curve d — u has been arbitrarily set to as described 
in the text. The errors are statistical only. 



the analysis of the full data set ||7I|. These data are in good qualitative agreement with the high-mass 
data published earlier |70|. The DY cross section ratio per nucleon for p + d to that for p + p is shown 
in Fig. 1^ as a function of X2. The acceptance of the spectrometer was largest for Xi;' = xi — X2 > 0. In 
this kinematic regime the DY cross section is dominated by the annihilation of a beam quark with a 
target antiquark. To a very good approximation the DY cross section ratio at positive xj? is given as 



(ynvip + d)/2aDY{p + p) ~ (1 + d{x2)/u{x2))/2. 



(12) 



In the case that d = u, the ratio is 1. Figure ^ shows that the DY cross section per nucleon for p + d 
clearly exceeds p + p, and it indicates an excess of d with respect to u over an appreciable range in X2- 
Figure |^ also compares the data with next-to-leading order (NLO) calculations of the cross section 
ratio using the CTEQ4M [^] and the MRS(R2) parton distributions. The data are in reasonable 
agreement with the predictions for X2 < 0.15. Above X2 = 0.15 the data lie well below the CTEQ4M 
and the MRS(R2) values. Following the publication of the E866 high-mass data several new 
parametrizations for parton distributions have been put forward. In Fig. ^ we compare the full E866 
data with the calculations using CTEQ5M |7|, MRST |6^, and GRV98 |75[. Also shown in Fig. | is 
a calculation using a modified CTEQ5M parton distributions. The modified CTEQ5M parton distri- 
butions, in which the d + u parametrization was maintained but d was set identical to u, were used to 
illustrate the cross section ratio expected for a symmetric d/u sea. The E866 data clearly show that 
d u. 
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Figure 6: The ratio of dju in the proton as a function of x extracted from the Fermilab E866 
cross section ratio. The curves are parametrizations of various parton distribution functions. 
The error bars indicate statistical errors only. Also shown is the result from NA51, plotted 
as an open box. 



Using an iterative procedure described in |7y, |7T], 0, values for dju were extracted by the E866 
collaboration at = 54 GeV^ and shown in Fig. At a; < 0.15, d/u increases linearly with x and is 
in good agreement with the CTEQ4M and MRS(R2) parametrization. However, a distinct feature of 
the data, not seen in either parametrization of the parton distributions, is the rapid decrease towards 
unity of the d/u ratio beyond X2 = 0.2. Figure ^ shows that the most recent parton distribution 
parametrizations (MRST, CTQ5M, GRV98) adequately describe the E866 data. The result from NA51 
is also shown in Figure ^ for comparison. 

The d/u ratios measured in E866, together with the CTEQ5M values for d + u, were used to obtain 
d~u over the region 0.02 < x < 0.345 (Fig. |^. As a flavor non-singlet quantity, d{x)—u{x) is decoupled 
from the effects of the gluons and is a direct measure of the contribution from non-perturbative processes, 
since perturbative processes cannot cause a significant d, u difference. From the results shown in Fig. |^, 
one can obtain an independent determination |7T| of the integral of d{x) — u{x) and compare it with 
the NMC result (Eq. g). E866 obtains a value /q [d{x) - u{x)] dx = 0.118 ± 0.012, which is 4/5 the 
value deduced by NMC |g. 

The E866 data also allow the first determination of the difference of the momentum frac- 
tion carried by d and u. One obtains Jo;o2^ x \d{x) - u{x)] dx = 0.0065 ± 0.0010 at = 54 GeV^. 
If CTEQ4M is used to estimate the contributions from the unmeasured x regions, one finds that 
Jq X d{x) — u{x) dx = 0.0075 ± 0.0011. Unlike the integral of d{x) — u{x), the momentum integral is 

dependent and decreases as increases. The dependence of the momentum fraction carried by 
various partons are shown in Fig. |^. The calculation uses both the MRS(R2) and the new MRST ^7 
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Figure 7: Comparison of the E866 d — u results at = 54 GeV^ with the parametriza- 
tions of various parton distribution functions. The data from HERMES |p2| are also shown. 



parton distributions for comparison. Figure ^ shows that the momentum fractions carried by the sea 
quarks increase with Q^. In contrast, the difference of the momentum fraction carried by up and down 
sea quarks decreases with Q^. 

3.5 HERMES Semi- Inclusive Experiment 

It has been recognized for some time that semi-inclusive DIS could be used to extract the flavor de- 



pendence of the valence quark distributions [77 . From quark-parton model, the semi-inclusive cross 



section, a];^, for producing a hadron on a nucleon is given by 

1 daUx,z) _ j:^eU,{x)Dl{z) 



j 6j fi {x) 



(13) 



where Df{z) is the fragmentation function signifying the probability for a quark of flavor i fragmenting 
into a hadron h carrying a fraction z of the initial quark momentum. Cj and fi are the charge and the 
distribution function of quark i, and (Jn{x) is the inclusive DIS cross section. Assuming charge symmetry 
for the fragmentation functions and parton distribution functions, one can derive the relationship 



d,,{x) _ 4i?^(x) + 1 
Uy{x) ~ A + R-^ix) ' 



(14) 



where 



K^{x) = {da'!^^ {x, z)/dz — da^^ {x, z)/dz)/{dap^ {x, z)/dz — da^ {x,z)/dz). 



(15) 
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Figure 8: Q-dependence of the proton's momentum carried by various partons calculated 
using the MRS(R2) and MRST parton distribution functions. 



Based on a large number of semi-inclusive charged-hadron events in muon DIS from hydrogen and 
deuterium targets, EMC extracted []78[ the values of dy{x)/uy{x) over the range 0.028 < x < 0.66. 
The EMC result agrees with neutrino measurements [[791 , |80| , and it demonstrates the usefulness of 
semi-inclusive measurements for extracting valence quark distributions. 

Soon after the report of GSR violation by the NMC, Levelt, Mulders and Schreiber pointed out 
that semi-inclusive DIS could also be used to study the flavor dependence of sea quarks. In particular. 



where 



d{x) — u{x) J{z)[l — r{x, z)] — [1 + r{x, z)] 
u{x) — d{x) J{z)[l — r(x, z)] + [1 + r{x, z)] 



(16) 



dap {x,z)/dz — da^ {x,z)/dz 
dall^ {x, z)/dz — da'^ {x, z)/dz' 



J{z) 



3 l + Dl 
5 1-AT" 



{^)/Dt{z) 



{z)/Dl^{zy 



(17) 



Unlike the situation for dy{x)/uv{x) which is completely independent of the fragmentation functions. 
Equation |TB| shows that fragmentation functions are needed to extract the values of d{x) — u{x). 



The HERMES collaboration [£2[ at DESY recently reported their measurements of charged hadrons 
produced in the scattering of a 27.5 GeV positron beam on internal hydrogen, deuterium, and ^He 
target. The fragmentation functions [z) were extracted from the '^He data, while the hydrogen and 
deuterium data allowed a determination of r{x, z). The values of {d — u)/{u — d) show no z dependence 
and are positive over the region 0.02 < x < 0.3, showing clearly an excess of d over u. Using the 
GRV94 LO parametrization of u{x) — d{x), the HERMES collaboration obtained d{x) — u{x) 



as shown in Fig. |^. The integral oi d — u over the measured x region gives /Q''g^2 '^(^ 



u[x\ 



dx 
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Figure 9: Comparison between the x{d — u) results from E866 with the parametrizations of 
various parton distribution functions. The dashed (sohd) curves correspond to PDFs before 
(after) the E866 results were obtained. 



0.107 ± O.Q21{stat) ± 0.017{syst). The total integral over all x is extrapolated to be 0.16 ± 0.03, 
consistent with the result from NMC It is gratifying that the results from E866 and HERMES are in 
rather good agreement, even though these two experiments use very different methods and cover very 
different values ((Q^) = 54 GeV^ in E866 and (Q^) = 2.3 GeV^ in HERMES). 

It should be mentioned that semi-inclusive DIS could be extended ^ to situations involving po- 
larized lepton beam and polarized targets in order to study the flavor dependence of the spin-dependent 
structure functions. Both the Spin Muon Collaboration (SMC) and the HERMES Collaboration |^ 
have reported the polarized valence quark distributions, Auy{x) and Adjj{x), and the non-strange sea- 
quark polarization, Aq{x). 



3.6 Impact on the Parton Distribution Functions 



After the evidence for a flavor asymmetric sea was reported by the NMC and NA51, several groups |72 



7^, |8^ performing global analysis of parton distribution functions all required d to be different from u. 



The NMC result constrained the integral of d — m to be 0.149 ± 0.039, while the NA51 result requires 
d/u to be 1.96 ± 0.13 at x = 0.18. Clearly, the x dependences of d — u and d/u were undetermined. 
Figures ^ and |^ compare the E866 measurements of d/u and x{d — u) with the parametrizations of the 
MRS(R2) and CTEQ4M [^, two of the most frequently used PDF's prior to E866's measurement. 

Recently, several PDF groups published new parametrizations taking into account of new exper- 
imental information including the E866 data. The parametrization of the x dependence of d — m is 
now strongly constrained by the E866 and HERMES data. In particular, d{x) — u{x) or d{x)/u{x) are 
parametrized as follows; 
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MRST 1^: 

d{x) - u{x) = 1.29x°-^^^(l - + 9.987X - 33.34x2), Ql = 1 GeV'^, 

GRV98 [^: 

d{x) - u{x) = 0.20x-°-"(l - x^^-^l - 13.3x°-^ + 60.0x), Ql = 0.4 GeV\ 
CTEQ5M 17|, H: 

d{x)/u{x) = 1 - 1.095X + 3.159x°-V[l + (x - 0.188)V0.01346], Ql = 1 GeV^. (18) 

As shown in Fig. these new parametrizations give significantly different shape for d/u at x > 0.15 
compared to previous parametrizations. Table || also lists the values of d — u integral from various 
various experiments and from recent PDF's. 

Table 1: Values of the integral J^{d{x) — u(x))dx from various experiments and parton distribution 
functions. 



Experiment /PDF 


Integral 


NMC 


0.148 ±0.039 


E866 


0.118 ±0.012 


HERMES 


0.16 ±0.03 


CTEQ4M 


0.108 


MRS(R2) 


0.162 


GRV94 


0.163 


CTEQ5M 


0.124 


MRST 


0.102 


GRV98 


0.126 



It is interesting to note that the E866 data also affect the parametrization of the valence-quark 
distributions. Figure |10| shows the NMC data for F2 — -F^ at = 4 GeV^, together with the fits of 
MRS(R2) and MRST. It is instructive to decompose -Ffl^) ~^2{^) i^^o contributions from valence and 
sea quarks: 



1 2 r - 

(x) — -F^(x) = -X [m„(x) — (i„(x)] ± -X u{x) — d{x) 



(19) 
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As shown in Fig. 

F2 — In order to preserve the fit to 



the E866 data provide a direct determination of the sea-quark contribution to 

F2, the MRST's parametrization for the valence-quark 



distributions, 



d^, is significantly lowered in the region x > 0.01. Indeed, one of the major new 



features of MRST is that d^ is now significantly higher than before at x > 0.01. Although the authors 
of MRST attribute this to the new 1^- asymmetry data from CDF |39| and the new NMC results on 
F2/F2 it appears that the new information on d{x) — u{x) has a direct impact on the valence-quark 
distributions too. 

Another implication of the E866 data is on the behavior of — F^ at small x. In order to satisfy 
the constraint /q^[m„(x) — dy{x)]dx = 1, the MRST values of m^(x) — dy{x) at x < 0.01 are now much 
larger than in MRS(R2), since Mt,(x) — dy{x) at x > 0.01 are smaller than before. As a consequence, 
Ff — F2 is increased at small x and MRST predicts a large contribution to the Gottfried integral from 
the small-x (x < 0.004) region, as shown in Fig. |ll|. If the MRST parametrization for Ff — F^ at 
X < 0.004 were used, NMC would have deduced a value of 0.252 for the Gottfried integral, which would 
imply a value of 0.122 for the d — u integral. This would bring excellent agreement between the E866 
and the NMC results on the d — u integral. 
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Figure 10: — as measured by NMC at = 4 GeV^ compared with predictions based 
on the MRS(R2) and MRST parametrizations. Also shown are the E866/NuSea results, 
rescaled to = 4 GeV^, for the sea-quark contribution to — -Pf . For each prediction, 
the top (bottom) curve is the valence (sea) contribution and the middle curve is the sum of 
the two. 



3.7 Comparison Between Various Measurements 

Are the various measurements of sea quark flavor asymmetry consistent among them? In particular, is 
the E866 result consistent with the earlier E772 and the NA51 DY experiments, and with the HERMES 
and NMC DIS measurements? To address this question, we first compare E866 with E772, both of which 
are DY experiments using 800 GeV proton beam with essentially the same spectrometer. Figs. |^ and ^ 
show that the MRST parton distributions, which determined d — u based on the E866 data, can also 
describe the E772 data very well, and we conclude that the E866 and E772 results are consistent. 

Although both NA51 and E866 measured a{p+d) /2a{p+p) to extract the values of d/u, some notable 
differences exist. As mentioned earlier, NA51 measured the ratio at a single value of X2 {x2 = 0.18) near 

~ using a 450 GeV proton beam, while E866 used an 800 GeV proton beam to cover a broader 
range of X2 at xp > 0. It is instructive to compare the NA51 result at X2 = 0.18 with the E866 data 
at X2 = 0.182. Table lists the kinematic variables and physics quantities derived from these two data 
points. It is interesting to note that the values of a{p + d)/2a{p + p) at X2 = 0.18 are actually very 
similar for NA51 and E866, even though the derived values for d/u differ significantly. This reflects 
the difference in for both experiments, making the values oi d/u extracted from a{j) + d) /2a{p + p) 
different. The other difference is Q^, being a factor of 3.6 higher for E866. Using MRST and CTEQ5M 
to estimate the dependence of d/u, we find that the NA51 value of d/u is reduced by ~ 3 % going 
from = 27.2 GeV^ to = 98.0 GeV^. This brings slightly better agreement bewteen NA51 and 
E866. 

The methods used by HERMES and E866 to determine d — u are different, and it is reassuring that 
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Figure 11: — -F2" measured by NMC at = 4 GeV^ compared with parametrization of 
MRST. The dashed curve corresponds to 0.21x^'^'^, a parametrization assumed by the NMC 
for the unmeasured small-x region when the Gottfried integral was evaluated. 



Table 2: Comparison between the 450 GeV NA51 result and the 800 GeV E866 data point |70 
X2 = 0.18. 



near 





(^2) 


{xf) 






d/u 


NA51 


0.18 


0.0 


5.2 


1.099 ± 0.039 


1.96 ± 0.246 


E866 


0.182 


0.192 


9.9 


1.091 ± 0.044 


1.41 ± 0.146 



the results came out alike, as shown in Fig. 0. The d~u values from HERMES are in general somewhat 
larger than those of E866. At a relatively low mean of 2.3 GeV^, the HERMES experiment could be 
subject to high-twist effects [Q. Additional data from HERMES are expected to improve the statistical 
accuracy. 

The comparison between E866 and NMC in terms of the integral oi d — u has been discussed earlier. 
A possible origin for the apparent differences of the integral was also discussed in Section 3.6. 

4 Origins of the d/u Asymmetry 

The earliest DIS experiments indicated that the Gottfried integral was less than 1/3, leading to specu- 
lation regarding the origin of this reduction. Field and Feynman suggested that it could be due to 
Pauli blocking in so far as uu pairs would be suppressed relative to dd pairs because of the presence of 
two u-quarks in proton as compared to a single rf-quark. Ross and Sachrajda questioned that this 
effect would be appreciable because of the large phase-space available to the created qq pairs. They 
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also showed that perturbative QCD would not produce a. d, u asymmetry. Steffens and Thomas |^ 
recently looked into this issue, explicitly examining the consequences of Pauli blocking. They similarly 
concluded that the blocking effects were small, particularly when the antiquark is in a virtual meson. 

The small d, u mass difference (actually, md > my) of 2 to 4 MeV compared to the nucleon 
confinement scale of 200 MeV does not permit any appreciable difference in their relative production 
by gluons. At any rate, one observes a surplus of d which is the heavier of the two species. As pointed 
out above, blocking effects arising from the Pauli exclusion principle should also have little effect. Thus 
another, presumably non-perturbative, mechanism must be found to account for the large measured 
J, u asymmetry. 

As many of the non-perturbative approaches to explain the asymmetry involve the use of isovector 
mesons (particularly pions), we present some of the motivation for believing that pion field is intrinsic 
to the nucleon's makeup. 

4.1 The Nucleon's Mesonic Field in Strong and Electroweak Process 

The requirement of the pion as intrinsic to the nucleon has been evident to most who have addressed 
the nucleon's interactions. To address the nucleon's partonic structure as totally independent of its 
interactions seems naive indeed. However, the success of the constituent quark model in simply and 
directly explicating nucleon properties has focused attention on the quark structure of the nucleon to 
the point that the role of pions was obscured. The magnetic moments of baryons are well characterized 
by constituent quarks, and pions are now known to have little effect on the result [^. In its most 



extreme form there was a hope to start with three constituent quarks at some low scale and simply 
employ QCD evolution to generate nucleon's observed partonic distributions. We now know this is 
most unlikely if not impossible, as non-singlet quantities such as qa, the u, d asymmetry, are of a 
non-perturbative origin and must be inserted ab initio if the nucleon's parton distributions are to be 
properly characterized. In the following we hst properties of the nucleon that require it to possess a 
meson cloud. 

• The nucleon's strong interactions, particularly the long-range part of the nucleon-nucleon inter- 
action have been characterized via meson exchange. The development of a low energy nucleon- 
nucleon potential has gone on for many years [^, ^ with the long-range part in particular 
requiring a dominant role for pion exchange. Attempts to generate this interaction from QCD- 



inspired models |9^, Q have not met with quantitative success so need for meson-exchange to 
account for the medium- and long-range parts of the nucleon-nucleon interaction appears beyond 
doubt. 

The requirement that the nucleon axial current be partially conserved (PCAC) requires the pion 
to be an active participant in the nucleon. As the pion is the axial charge, it has a dominant status 
in PCAC. Figure [T^ is a diagram of the pion's role in nucleon beta decay. Employing PCAC, the 
Goldberger-Treiman relation can easily be derived, 

9A = (20) 

In the above expression F^^ is the pion decay constant (92.42 ± 0.26 MeV), Qc is the charged 
pion (iinp) coupling constant (47r(14.17 ± 0.2)^/^) and Mp the proton mass. This yields a 



value for gA that is (3.8 ± 2.5) % too high, not inconsistent with the variance expected due 
to the breaking of chiral symmetry {mu-,rnd-,rnTj ^ 0). Furthermore, the value of the induced 
pseudoscalar form factor, QpiQ"^), is also directly dependent on the pionic field of the nucleon. 
PCAC yields a result [|og] 

9p{Ql = -0.88mJ) = 8.44 ± 0.23, (21) 
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Figure 12: Diagram showing pion's role in nucleon beta decay. 



consistent, albeit with large errors, with the measured value, 

(7p(Q2) = 8.7±2.9. (22) 

There are of course many other examples where mesons, especially the pion are employed to 
account for the nucleon's properties. Indeed, QCD based chiral models have treated the pion as 



a basic degree of freedom from the earliest times |101|. However, the consequences of their effects 



were not incorporated into parton distributions until the early 90's [|102| , |103|| . 



As is well known, the information on the nucleon's structure from parity-conserving elastic electron 
scattering is contained in four form factors. They are the Dirac (Fi(Q^)) and Pauli {F2{Q'^)) form 
factors for the neutron and proton. Linear combinations of these form factors are often employed, for 
example an isoscalar [F^ + F") and an isovector [F^ — F") combination as well as the Sachs | 1U4 
electric and magnetic form factors: 

= - ^,F,{Q')- Gm{Q') = F,{Q') + F,{Q'). (23) 

The normalization of the form factors are Ff (0) = 1, F^{Q) = 0, Fi{0) = 1.793, and F2''(0) = -1.913. 
The Sachs form factors are just the Fourier transforms of the charge and magnetization distributions 
in the Breit frame. 

The large body of electron scattering data spanning < < 30 GeV^ can be analyzed in a largely 
model independent fashion using dispersion relations. They fit the data rather well but are not very 
revealing of the underlying structure of the nucleon. The electric form factor of the neutron on the 
other hand could be most revealing regarding the presence of a pvr" component in the neutron wave 
function. There are two difficulties to be faced. First, the presence of the Foldy term in G^(Q^) tends 
to obscure and complicate the interpretation of the neutron charge distribution ||105|| . Secondly, there 



are QCD hyperfine interactions ||105| , |106| , [L07|| that produce a negative neutron charge radius similar 
to the p7T~ component. Thus it will take extensive additional theoretical and experimental work to 
clarify this important issue in electron-nucleon scattering. However, some promising steps are being 
taken [M . 

An area where appreciable recent progress has been made is in the analysis of the nucleon's electro- 



magnetic form factors in terms of relativistic baryon chiral perturbation theory | 109 |. This approach is 
the effective theory of QCD at low energy, which limits its range of applicability to ~ < O.SGeV^. 
The vector mesons are brought in without the introduction of additional parameters by using the results 
of theoretical dispersion analysis of the form factors. The neutron charge form factor is little changed 
by the inclusion of vector mesons, while they have a beneficial effect on the three other nucleon form 
factors. Figure ^ shows the level of agreement obtained in G^(Q^) applying relativistic baryon chiral 



19 



0.15 




Figure 13: Comparison between data and calculations ||109|| using relativistic baryon chiral 
perturbation theory for G%{Q'^). The dashed curves are the results of calculation without 
vector mesons, the upper dashed curve being to 3rd order while the lower is to 4th order. 
The solid curves are the results including vector mesons where the upper curve is to 3rd 
order while the lower is to 4th order. The dot-dashed curve is the result of the dispersion 
theoretical analysis. 



perturbation theory as shown in Ref. 109 . While the fits are impressive it still remains illusive as to 
what is added to our picture of the nucleon. The pion field determines most of the properties of the 
form factors. For example, the leading terms for the radii of the isovector charge and magnetic form 
factors ||108|| are 



where is the nucleon anomalous isovector magnetic moment (3.706). Note that in the chiral limit 
(Mtt — > 0), both radii become infinite, a not surprising result for massless pions. Our main purpose 
is to show that the pion mass is a critical parameter in the nucleon form factors in this chiral model 
derived from QCD. The important role of pions in the structure of the nucleon appears to us to be 
beyond question and the nucleon's parton distributions should reflect this. 

4.2 Meson-Cloud Models 

As the u, d asymmetry cannot be generated via perturbative process, several non-perturbative models 
have been proposed that yield an asymmetry. The most straightforward of these models are those 
that attribute the asymmetry to the existence of a "pion cloud" in the proton. The relevance of pion 
cloud for sea-quark distributions appears to have first been made by Thomas in a publication ||110 | 



treating SU(3) symmetry breaking in the nucleon sea. Sullivan | 111 | had earlier shown that virtual 



meson-baryon states directly contribute to the nucleon's structure function. A large number of authors 



have contributed to calculating the asymmetry from this perspective, so recent reviews ||112| , |113|| should 
be consulted for a complete list of contributions. 

In the meson-cloud model, the virtual pion is emitted by the proton and the intermediate state is 
pion + baryon. More specifically, the proton is taken to a linear combination of a "bare" proton plus 
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pion-nucleon and pion-delta states, as below, 



VI - a - b\po) + M-\ll\PoA + /f 1%^^)) + ^K\ll\^>-) - + ^|A>+)).(25) 



\P) 

The subscript zeros on the virtual baryon states indicate that they are assumed to have symmetric 
seas, so the asymmetry in the antiquarks must be generated from the pion valence distribution. The 
coefficients a and b are the fractions of the ttN and vrA configurations, respectively, in the proton. 
These fractions can be calculated using the irNN and vriVA couplings, and form factors as taken from 
experiment. The asymmetry in the proton sea arises because of the dominance of vr"*" among the virtual 
configurations. These calculations, to be discussed in detail below, reproduce many aspects of the data 
but suffer two problems. First, there is a strong dependence on the value used to cut off the integral over 
the form factor, and secondly, any such calculation of d{x)/u{x) is unreliable because the magnitude of 
the symmetric sea is unknown. Before launching into a more detailed presentation of these calculation, 
a few general observations can be made. It is instructive to examine some general properties of the irN, 
ttA ansatz. A useful expression to consider is 

dpjx) ^ |a/i^(x) + |&/,^(3:) + ig(a;) 

u,{x) laf^{x) + lbf^{x) + lS{xy ^ ' 

S{x) is the amount of symmetric sea and f^{x){f^{x)) are functions that characterize the distributions 
in X associated with the irN (vrA) configurations. From this expression it is clear that the magnitude 
of the symmetric sea must be known if the ratio is to be predicted. The maximum value that the ratio 
can assume is 5, if there were no symmetric sea and only vrA^ configurations contribute. The minimum 
value of the ratio is 1/2 which occurs in the absence of a symmetric sea and only vrA configurations 
contributing. The ratio takes on a value 1 with a pure symmetric sea, or a symmetric sea with a = b/2 
and fj^ = f^. Pion models are much more effective in dealing with the integral isolating the contribution 
from asymmetric sea (AS) 

/■I - _ 1 

Ias = / [dp{x) - Up{x)]dx = -(2a - b), (27) 
Jo o 

as there is no contribution from the symmetric sea in this case. 

Many attempts |p.l4| , |115| , |116| , |117| , |118| , |119| , |120|| have been made to calculate the flavor asymmetry 



due to isovector mesons. Most start with the following convolution expressions: 

xqp{x,Q'^) = X1«MB / dy fuBiv) - qni-^Q^), (28) 

fr^ Jx y y 



MB 



where 



fMsiy) = -T^y / dt— 2-^^Mpb(^> A), t^in = rripV - m^- . (29) 

IDTT J —CO \t — ) — y 

In the above expressions x is the fraction of proton's momentum carried by the antiquark, and y is 
fraction carried by the meson (M). The meson-proton-baryon couplings are characterized by coupling 
constants gupBi and form factors FMpsit, A) where A is a cutoff parameter. F{t, rrip, rriB) is a kinematic 
factor depending on whether B is in the baryon octet or decuplet. As pions are the only mesons usually 
considered and the baryons are usually restricted to nucleons and deltas, the coupling constants are 
well known and the partonic structure of the pion, qn{x, Q"^), is flxed by measurement of the DY process 
using high energy pion beams. The only uncertainties are the form factors F^pAr(t) and F^pa(^)- One 



attempts to determine these form factors by using ||121|, |122|| the measured yields from a variety of high 
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Figure 14: Comparison of the E866 d—u results at = 54 GeV^/c^ with the predictions 
of pion-cloud and chiral models as described in the text. 



energy hadronic reactions at small pt such as pp — > nX and pp Even though there is a 

sizable amount of available data, employing such a procedure does not produce a precise result. The 
cutoff parameters used in the extracted form factors are of the order of 1 GeV but the uncertainties in 
their values produce factors of two in the predicted antiquark content of the nucleon due to pions. 

Even though calculated value of the integral of d{x) — u{x) is often in agreement with experiment, 
it is more difficult to achieve a quantitative fit to the measured x dependence of the difference . 

■u{x) from E866 with a pion-cloud model calculation, following the procedure 

1.0 GeV for the vrA^A^ form factor and A = 0.8 



Figure [Tj compares d{x 

detailed by Kumano ||115|| . A dipole form, with A 
GeV for the vrA^A form factor, was used. Calculations of d{x)/u{x) are even more unsuccessful, as 
knowledge of the x dependence of the symmetric sea is required in this instance [j76|. 

It is instructive to compare the pion-model prediction with the current PDF parametrization of 
x{d + u). Figure |15| shows that at small x {x < 0.1) the valence quarks in the pion cloud account for 
less than 1/3 of the d + u content in the proton. In contrast, at large x {x > 0.5) the pion model would 
attribute all of d + u to the pion cloud. 

Determining the appropriate cutoff is difficult because the high energy hadronic reactions {pp —>■ irX, 
pp — >■ nX, pp — ^ A^+X) used to fix the form factors all have experimental backgrounds, and extracting 
quantitative results from such reactions is difficult. Nikolaev et al. ||122|| carefully investigated these 
reactions by including constraints from the Regge behavior of various mesons on the total photo- 
absorption cross section. They find that the contribution of p and 02 Reggeons to the proton structure 
function is negligble. They also find the relative contribution of vrA to be much smaller than previous 
analysis. For example a = 0.105,6 = 0.015 for cutoffs of Rq = 1.5GeV~^ and Rq = 2 GeV~^ for the 
ttN and vrA respectively. Their fits to dp( 
In a recent publication 



are shown in Fig. 16. 



123|] Thomas et al. appear to have avoided the form factor issue by having 
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Figure 15: Comaprison of x{d + u) obtained from the valence quarks in the pion cloud with 
the parametrization of the MRST parton distribution functions. 



recourse to Chiral Perturbation Theory [ 124|| . In this case the theory is implemented with pions, 



nucleons and deltas as the degrees of freedom. They find for the O'th moment of the dp 
to the vrA^ loop to be 



I^^\n= f dx[d{x) - u{x)] 

J 



-mln 



ml 



due 



(30) 



where qa is the axial charge of the nucleon, the pion decay constant and /i a mass parameter. They 
stress that this leading non analytic (LNA) behavior originates from one-pion loop and hence is robust. 
General expressions are derived for the n'th order terms which are suppressed by powers of (m^/Mp)". 



Using Eqs. (2), (5), (6) and (8) from Ref. ||123|| one obtains the ratio of the contribution of ttA to that 
of TiN to be 



tLNA 






ttA 


tLNA 




^AS 


■kN 



3 (Ma + 



N) 



Evaluating the net result using /z 



75 

47r/^ one finds 



Ml 



-0.12. 



J^^^ = 0.154 



(31) 



(32) 



close to the E866 value of 0.118±0.012. 

The values for I as that result from the various parametrizations of the nN and ttA virtuality are 
all reasonable and not far from the experimentally measured value of 0.118ib0.012. The values for a 
and h from various authors who have employed the meson cloud model show that a typical value of 
h/a is 1/2. This leads to a = 0.24, h = 0.12 to satisfy the observed flavor asymmetry, resulting in the 
probability of finding a pion in a nucleon (Eq. |2^) of a + 6 = 0.36. 
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The ttN , 71 A model also allows one to evaluate the effect of virtual meson emission on the spin 
carried by quarks. The pion is emitted from the proton in a p-wave and therefore reduces the spin 
projection of the virtual baryon, thus reducing the spin projection carried by quarks. In this picture 
there is no spin on the antiquarks and the reduced spin in the baryon sector is carried off in the orbital 
motion of the pions. The expressions below for the spin on the quarks assumes an SU(6) wave function 
for the baryons. 



Au + Ad=l 



2, 

3(2a 



Has 



(33) 



Where 



. , 5 20, ,,32 

Au- Ad= (2a + 6) H 

3 2T ' 21 



Aq = [ \q{x) t +q{x) t -qix) | -q{x) l]dx. 
Jo 



(34) 



(35) 



Interference occurs in Eq. ^ because the operator aztz connects nucleon and delta states, while 
in Eq. |3^ the operator az cannot. Thus, virtual pion emission is seen to reduce the spin carried by 
quarks. With no further assumptions we have Am + Ad = 0.76 ± 0.024, which is appreciably smaller 
than one but still greater than the observed value of 0.40 ± 0.02. Employing a commonly used value 
of 6 ~ a/2, we obtain Au — Ad = 1.51 ± 0.053. Again virtual pion emission reduces Am — Ad = qa 
from its SU(6) value of 5/3 but not nearly to the measured value of 1.256. Interestingly, the value 
extracted for qa using the measured flavor asymmetry in this model is nearly identical to that obtained 
by Weinberg |125| using chiral perturbtion theory and current algebra techniques. 



4.3 Chiral Models 



An alternative approach [||, |T2|, |123, [T2|, [T2|, [T30|, [131 



the d^u asymmetry is usually referred to as chiral models jlOl 



also employing virtual pions to produce 
After the pubhcation of the NMC 
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result on the violation of the Gottfried sum rule, Eichten, Hinchliffe and Quigg (EHQ) following 
an earlier suggestion of Bjorken [|132|| , published an interesting approach to explain the NMC result. 



They used chiral perturbation theory, employing the Manohar and Georgi ||101|| model in which the 



degrees of freedom are Goldstone Bosons (GB), and constituent quarks (i.e., U and D for up and down 
constituent quarks). The pion in this approach is emitted by a constituent quark. In this model most 
of the parameters are more or less prescribed (the mass of the constituent quark must be assigned). If 
a is the probability of a f/ constituent quark becoming tt^D, then to first order in a, invoking isospin 
invariance and the conservation of probability 

U {l-^a)U + a7T+D + ^an^U (36) 

with a corresponding expression for the down constituent quark. As the proton is UUD, to first order 
in a one has 

7(x 11a 

p^2U + D + —{u + u) + ^{d + d). (37) 

Thus the integrated asymmetry, I^b, is equal to a. Note that in the absence of a symmetric sea the 
maximum value for d{x)/u{x) in this model is 11/7. Restricting themselves to only the pion the authors 
calculate a value of a = 0.083. This value is somewhat small, but within 50% of the measured value. 
However, it is very interesting that their calculation of a, carried out in the context of chiral field 
theory ||101|| generates a probability of finding a pion in a nucleon of 9a/2 = 0.374, very similar to that 



found in the meson-cloud model. 

The authors also investigated the consequences of this model on the quark spins. Using SU(6) wave 
function for the proton and assuming total spin-flip of the constituent quark upon emitting a pion, they 
found for the modified quark spins in the proton 

Au + Ad = l-3a = l- 3Ias (38) 

Au-Ad = -- -a. (39) 
3 3 

Using a = 0.118±0.012 from the E866 measured flavor asymmetry, Au+Ad = 0.65 and Au—Ad = 1.47. 
Both are reduced from the SU(6) values but are not as small as the measured values of 0.40 and 1.256, 
respectively. One believes that the SU(6) value is too large due to relativistic effects but the size of 
reduction is highly model dependent. 

The differences that occur between the EHQ model and the meson-cloud model employed above are 
that the pion is emitted by a single constituent quark rather than by the proton. EHQ have a single 
parameter a which corresponds to fixing a/b at 5/4 in the meson-cloud model. We believe that this 
value is too small and does not easily produce the large values of d{x)/u{x) that are observed in the 
interval 0.1 < x < 0.2. EHQ achieves greater suppression of the quark spin because they assume a total 
spin flip of the constituent quark upon emission of a pion. 

Reference |^6|] also examined the effect of extending the space of GBs beyond the pion. Treating the 
octet (vr, K, rf) of GBs as degenerate with equal coupling to all members, the proton becomes to first 
order in a', 

p^2U + D + 2a'{u + u) + -a'{d + d) + —a'{s + s). (40) 

3 3 

In this case the u, d flavor asymmetry is diluted, with the maximum value for d{x)/u{x) becoming 1.333, 
and Ias = 2a' /3. Thus satisfying the E866 measured value for Ias requires a' = 0.177 ± 0.017. The 
effect on the proton spin is given by 

A , 5 45a' 

Au- Ad= (41) 

3 9 ^ ' 
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Au + Ad + As = l 



16a' 



(42) 



Both larger value for a' and the larger number of mesons available to flip the quark spins conspire to 
produce a large effect on the spin values, Am — Ad = 0.978 ± 0.066 and Am + Ad + As = 0.056 ± 0.091. 
It is clear that this approach produces too little flavor asymmetry and too much reduction of quark 
spin. In addition one can see from Eq. ^ that far too much mesonic sea has to be generated to produce 
the observed flavor asymmetry. Continuing in this vein, EHQ showed that the u, d asymmetry vanishes 
in the case of an ideal U(3) nonet. This of course is not what is observed as the t] and r]' contributions 
are suppressed because they are much more massive than the vr. 

This approach has been considerably extended by Cheng and Li |p.26| , |127| , |133| , |134|| They employ 
the full meson U(3) nonet but allow the symmetry to be broken U{3) — >■ SU{3) x f/(l). The coupling 
to the singlet rj' is unique. Originally they made an ad hoc assumption that the singlet coupling was 
the negative of the octet coupling. Subsequently they have presented strong support for this point of 
view using the coupling of quarks to instantons ||135| , |136|| and using the same physics that produces the 
large mass for the rj'. Defining the ratio of the coupling to the singlet to that of the octet as fi/ fs = 
they found ( = —2. The ratio d/u in terms of (, due to the full nonet of GBs is 



d 



u 



3^ 



2 + 



(43) 



which gives 1 for = 1, and 4/3 for Q = (the pure octet case). With ( = —2 the ratio is 2, 
an altogether reasonable value. However, in their most recent tabulation [p.34|| of results they choose 
( = —1 and generate the results shown in the third column of Table ^. In order to account for expected 
suppression factors arising from mass differences, they included a weighting factor, {{k"^) +^gb)~^j 
each GB amplitude. This result is shown in column 4 and is seen to reduce the excessive contribution 
of strange quarks. This description while apparently quite successful, has the problem of producing a 
softer X distribution for d{x) — u{x) than is observed in E866. Figure |I^ shows a chiral-quark model 
calculation using the formulation described in Ref. ||129|| . The predicted d{x) — u{x) distribution 
is too soft because the GB is emitted by a constituent quark, which carries only 1/3 of the nucleon's 
momentum. This suggests that additional dynamics must be included in the chiral model if it is to 
produce agreement with the data. 



Table 3: Comparison of the measured parton flavor and spin structure of the proton with the model 
of Cheng and Li [|133| , |134|| . Some of the measured values are slightly different from the current values. 
The value of (fc^) used is 350 MeV^. 



Quantity 


Measured Value 


SU3 


Broken SU3 


d — u 


0.147 ±0.026 


0.15 


0.15 


2s/{u + d) 


~ 0.5 


1.86 


0.6 


Au 


0.82 ±0.02 


0.78 


0.85 


Ad 


-0.43 ±0.02 


-0.33 


-0.40 


As 


-0.10 ±0.02 


-0.11 


-0.07 


AS 


0.29 ±0.06 


0.34 


0.38 


Au,Ad 


0.01 ±0.07 


0.0 


0.0 


9A 


1.257 ±0.03 


1.12 


1.25 



By way of further comment on Table |^, there is no reason to expect the ratio of 2s/ {d + u) generated 
from GBs alone should be 0.5, as this ratio depends sensitively on the contribution from the symmetric 
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sea. Also there is a lack of consistency in introducing a mass factor in the GB sector and ignoring the 
effect in the baryon sector. The baryons accompanying the GBs are A^, A,A, S with mass splittings 
comparable to those occuring among the GB. 

Continuing with approaches that explain the flavor asymmetry and spin quenching by introducing 
the role of mesons, a novel approach has been taken by Ball et al. ||5^, ^ in which the mesons are intro- 
duced by including them in the evolution of the quark structure function functions as a non-singlet 
non-perturbative component. The approach is described in detail in pP[ and incorporates many 
attractive features. For example they explicitly account for the mass dependence of the contribution 
from each of the pseudoscalar mesons. Their approach generates a very strong scale dependence to 
the value observed for the Gottfried Sum Rule. Unfortunately, there is no relevant DIS data to verify 
their prediction. Drell-Yan experiments which could equivalently measure I as as a function of are 
typically restricted to using lepton pairs above the {Mq,i = 3.69 GeV), so have minimum of ~ 16 
GeV^ making it impossible to observe the predicted scale dependence. 



A very nice paper recently published by Henley, Renk, and Weise ||137|| showed that the distribution 
in X of d{x) — u{x) measured by E866 corresponds to that associated with a pion cloud. Their result 
is somehow implicit in earlier work that showed the pion cloud model reproduced the x distribution 
observed in experiment, however this work makes explicit that the spatial distribution involved is what 
one would expect for a pion cloud. The method employed uses the coordinate space representation 



developed earlier by Filler et al. | |138| | and Vanttinen et al. |139| . They introduce a light-cone dimen- 



sionless space-time variable z = y ■ P, where P is the nucleon momentum. The light-cone distance is 
= t + = 2z/M, with M the nucleon mass. The dimensionless variable is conjugate to Bjorken x, 
and z = 10 corresponds to y~^ = 4 fm. The authors work with parton distribution functions such as 
CTEQ5 or MRST that accurately capture the measured d{x) —u{x) distributions. The transformation 
from the x distribution is carried out via 

{D - U){z,Q^) = f\d{x,Q^) -u{x,Q^)]sin{zx)dx (44) 







Figure ^ shows how the (Z) — U){z) distributions appear when plotted against z. The peak in z is at 
2.5 or about Ifm and the bulk of the distribution is between z = 1 and 10 corresponding to 0.4 to 4 
fm, in line with what one expects for a p-wave meson cloud. 
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4.4 Other Approaches not Directly Including Mesons 



4.4.1 Instanton Models 



Instantons have been known as theoretical constructs since the seventies ||140| , |141| , |142|| . They represent 
non-perturbative fluctuations of the gauge fields that induce transitions between degenerate ground 
states of different topology. In the case of QCD, the collision between a quark and an instanton flips 
the helicity of the quark while creating a qq pair of different flavor. Thus, interaction between a u 
quark and an instanton results in a m quark of opposite helicity and either a dd or ss pair. Such a 
model has the possibility of accounting for both the flavor asjTiimetry and the "spin crisis" 
However, the model has proven difficult to exploit for this purpose. There is only one case 



143, 144 



T35I of its 



being directly employed to explain these anomalous effects. In the case of the c?, u flavor asymmetry, 
the authors of ref. [ 145|| fit the instanton parameters to reproduce the violation of the GSR observed by 
NMC The prediction ||145|| at large x, d{x)/u{x) —>■ 4, is grossly violated by experiment (see Fig. |^). 
Thus, it appears that while instantons have the possibility for accounting for flavor and spin anomalies, 
the approach is not yet sufficiently developed for a direct comparison. The final state created via an 
instanton collision is quite similar to that created via the emission of meson in the chiral model. 

However, it must be pointed out that the instanton model contains elements that strongly influence 
other descriptions. For example, as cited above in the work of Cheng et al. [ |135| , |136|] , instantons provide 
the basis for generating a flavor asymmetry when the full GB nonet is employed. If the coupling to the 
f]' is the same as the coupling to the GB octet, the sea is flavor symmetric. They originally made the 
ad hoc assumption that the rj' coupling was equal in magnitude but opposite in sign to the coupling 
to the octet. This assumption has been justified | |13(j| ] by invoking instantons. Further, as will be 
discussed below, the chiral quark soliton model which is developed using a vacuum sea of pions can 
equally well | |146[ be derived using the instanton model of the QCD vacuum. Thus it appears that the 
notion of instantons might be essential to the flavor asymmetry and that they play a crucial role in 
underpinning what ever explanation is used. 



4.4.2 Lattice Gauge Approach 

It would be extremely informative to explain the d, u asymmetry in terms of lattice gauge calculations 
as it would provide some insight as how the evolution of the asymmetry takes place. Indeed, soon after 
the existance of the asymmetry was established K.F. Liu and his collaborators investigated ||147| , |148 
the issue. Employing a path-integral formalism they established that the d, u difference comes from the 
set of diagrams they termed connected diagrams (Fig. ^ (a,b)). The other sources of sea quarks are 
disconnected diagrams, shown in Fig. |18](c). The loops in the disconnected diagrams are generated by 
gluons, and can produce all flavors of sea quarks. These diagrams produce equal numbers of up and 
down sea quarks. The loops in the connected diagrams can only produce up and down sea quarks and 
presumably are the source of the up, down asymmetry. 

This approach has the interesting feature that the strange quark sea is created entirely from discon- 
nected diagrams and hence the s and s distributions are identical, in good agreement with the results 
from neutrino DIS. It would also account for the near vanishing of the strange vector form factors for 
the nucleon as extracted from parity violating electron scattering. However it leaves unspecified the 
spin carried by strange quarks. 



4.4.3 Chiral-Quark Soliton Model 

One of the interesting approaches to emerge over the past 5 years is that of the chiral-quark soliton 
model ||146| , |149| , |151| , |152|| . It uses the large limit of QCD, which becomes an effective theory of 
mesons with the baryons appearing as solitions. At low energy the effective dynamics is described by a 
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(a) (b) (c) 



Figure 18: Euclidean path integral from Ref. ||148|| for evaluating the hadronic tensor W, 



UfJ,- 



(a) and (b) are connected insertions and (c) is a disconnected insertion. 



chiral Lagrangian for the pion, which appears as the GB of the spontaneously broken chiral symmetry. 
The model employs a realistic effective interaction, containing all orders of derivatives of the pion field, 
defined by the integral over quark fields having dynamically generated mass and interacting with the 
pion field in a minimally chirally invariant way. It is valid over a range of momenta up to the inverse 
of the instanton size = 0.6 GeV). The parton distribution obtained at this low momentum scale 
can then be evolved to higher for comparison to experiment. 

Quarks are described by single particle wave functions which are solutions of the Dirac equation 
in the field of the background pions. The spectrum of single particle states includes a single bound 
state level plus distorted positive and negative Dirac continumn. The discrete level and the negative 
continuum are occupied producing a state of unit baryon number. The distribution Au{x) — Ad{x) 
appears in leading order (A*"^) in a l/N^ expansion while u{x) — d{x) appears in next-to-leading order. 
For example the isovector spin carried by quarks is 

Au{x) - Ad{x) = --{2n)N,MN E J T^K{^){l + 'lVh5r'5{h + En + xMn)^4k). (45) 



n(occup) ' 



Where ^n(k) are the single particle wave functions, and 2T3 = ±1 for the proton and neutron respec- 
tively. The sum runs over all occupied single particle quark states, the bound state and the negative 



continuum. Figure |19| shows some of the results from Ref. 



The values for the integral of the quantities shown in Fig. [T^ are 



f [d{x) - u{x)]dx = 0.17; f [Ad{x) - Au{x)]dx = -0.31. (46) 

Jo 



Equation ^ presents a value for the integrated flavor asymmetry that is larger, but in the same ballpark 
as the measured value (0.118 ± 0.012). The distribution in x is not unlike that of models employing 
virtual mesons and hence should have a similar spatial distribution. However, the magnitude of the 
integrated isovector longitudinal-spin distribution for the antiquarks is surprisingly large and totally at 
variance with the results of models employing GBs to account for the flavor asymmetry. The GB models 
obviously require Ad{x) = Au{x) = 0. While one might think that the chiral soliton model should 
be consistent with the GB models as the baryonic soliton is developed in a background field of pions, 
the result for the spin carried by the antiquarks is totally at odds with any GB description. Hence 



one cannot view this non-pertubative approach as generating pseudoscalar mesons. Recent data |^ 
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Figure 19: Calculations from a cliiral-quark soliton model [p.53|| (a) The value calculated 
for d{x) — u{x) in the proton compared to early data from E866. The calculation has been 
evolved from Q = 0.6 GeV to Q = 7.35 GeV. (b) The calculated flavor asymmetry of the 
longitudinally polarized antiquark distribution in the proton (solid line) and the transversity 
distribution of the antiquarks (dashed line). Both distributions are aX Q = 0.6 GeV. 



showed Aq = 0.01 ± 0.04 ± 0.03, consistent with a small sea polarization. We await the judgement of 
future experiments to decide which description best conforms to reality. 



5 Further Implications of the Meson-Cloud Models 



5.1 Strange Sea of the Nucleon 

Models in which virtual mesons are admitted as degrees of freedom have implications that extend 
beyond the d, u flavor asymmetry addressed above. They create hidden strangeness in the nucleon via 
such virtual processes as p — *• A + K^, + K, etc. Such processes are of considerable interest as they 
imply different s and s parton distributions in the nucleon, a feature not found in gluonic production of 
ss pairs. This subject has received a fair amount of attention in the literature [ 121 , 1^ , 155 , |156| , |157|| 
but experiments have yet to clearly identify such a difference. Thus in contrast to the d, u flavor 
asymmetry, to date there is no positive experimental evidence for ss contributions to the nucleon from 



virtual meson-baryon states |19, 158 



A difference between the s and s distribution can be made manifest by direct measurement of the 
s and s parton distribution functions in DIS neutrino scattering, or in the measurement of the 
dependence of the strange quark contribution (-Ff^(q'^)) to the proton charge form factor. This latter 



case is not well known and follows from a suggestion of Kaplan and Manohar ||159|| regarding the new 
information contained in the weak neutral current form factors of the nucleon. Measurement of these 
form factors allows extraction of the strangeness contribution to the nucleon's charge and magnetic 
moment and axial form factors. The portion of the charge form factor Fig{q^) due to strangeness 
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clearly is zero at = 0, but if the s and s distributions are different the form factor becomes non-zero 
at finite q^. These "strange" form factors can be measured in neutrino elastic scattering ||160|| from the 



nucleon, or by selecting the parity- violating component of electron-nucleon elastic scattering, as is now 



being done at the Bates ||161|| and Jefferson Laboratories ||162 |. 

It is worth pointing out that there is a relationship between the parton distributions and the form 
factors of a hadron. If the neutron's charge form factor is explained in terms of a particular meson- 
baryon expansion, then one should expect that the expansion is consistant with the neutron's partonic 
structure. Little work appears to have been done bringing these descriptions together. Below is an 
example showing the impact of the nucleons pionic content, inferred from the flavor asymmetry, on the 
spin distribution in the nucleon. 

5.2 Sea Quark Distributions in Hyperons 

Dilepton production using meson or hyperon beams offers a means of determining parton distributions 
of these unstable hadrons. Many important features of nucleon parton distributions, such as the flavor 
structure and the nature of the non-perturbative sea, find their counterparts in mesons and hyperons. 
Information about meson and hyperon parton structure could provide valuable new insight into nucleon 
parton distributions. Furthermore, certain aspects of the nucleon structure, such as the strange quark 
content of the nucleon, could be probed with kaon beams. 

No data exist for hyperon-induced dilepton production. The observation of a large d/u asymmetry 
in the proton has motivated Alberg et al. ||163| , |164|| to consider the sea-quark distributions in the S. The 



meson-cloud model implies a d/u asymmetry in the S"*" even larger than that of the proton. However, 
the opposite effect is expected from SU(3) symmetry. Although relatively intense S"^ beams have been 
produced for recent experiments at Fermilab, this experiment appears to be very challenging because 
of large pion, kaon, and proton contaminations in the beam. 

6 Future Prospects 

6.1 d/u at Large and Small x 

The interplay between the perturbative and non-perturbative components of the nucleon sea remains to 
be better determined. Since the perturbative process gives a symmetric d/u while a non-perturbative 
process is needed to generate an asymmetric d/u sea, the relative importance of these two components 
is directly reflected in the d/u ratios. Thus, it would be very important to extend the DY measurements 
to kinematic regimes beyond the current limits. 

The new 120 GeV Fermilab Main Injector (FMI) and the proposed 50 GeV Japanese Hadron Facil- 
ity ||165|| (JHF) present opportunities for extending the d/u measurement to larger x {x > 0.25). For 



given values of xi and X2 the DY cross section is proportional to 1/s, hence the DY cross section at 



50 GeV is roughly 16 times greater than that at 800 GeV! Figure |20| shows the expected statistical 



accuracy for a{j) + d)/2a{p + p) at JHF ||166| , |167|| compared with the data from E866 and a proposed 



measurement [ p.68|| using the 120 GeV proton beam at the FMI. A definitive measurement of the d/u 
over the region 0.25 < a; < 0.7 could indeed be obtained at FMI and JHF. 

At the other end of the energy scale, RHIC will operate soon in the range 50 < ^/s < 500 
GeV/nucleon. The capability of accelerating and colliding a variety of beams from p + p, p + A, 
to A + A at RHIC offers a unique opportunity to extend the DY d/u measurement to very small x. 
Such information is important for an accurate determination of the integral of d — u, as well as for 
a better understanding of the origins for flavor asymmetry. The statistical accuracy for measuring 



a{p + d)/2a{p + p) in a two-month PHENIX run is shown in Fig. |2l|. Also shown in Fig. |2l| are the 
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Figure 20: Projected statistical accuracy for a{p + d)/2a{p+p) in a 100-day run at JHF [|166 



167|| . The E866 data and the projected sensitivity for a proposed measurement ||168|| at the 
120 GeV Fermilab Main-Injector are also shown. 



data from E866. 
in E866. 



The lowest X2 reachable at RHIC is around 10 ^, an order of magnitude lower than 



6.2 W Production 



To disentangle the d/u asymmetry from the possible charge-symmetry violation effect [|169|, |170|, |171 



one could consider W boson production, a generalized DY process, in p + p collision at RHIC. An 
interesting quantity to be measured is the ratio of the p + p ^ + X and p + p ^ W~ + X cross 
sections ||172|| . It can be shown that this ratio is very sensitive to d/u. An important feature of the 
W production asymmetry 'm p + p collision is that it is completely free from the assumption of charge 
symmetry. Figure ^ shows the predictions for p + p collision at -^i = 500 GeV. The dashed curve 
corresponds to the d/ u symmetric MRS SO' ||173|] structure functions, while the solid and dotted curves 
are for the d/u asymmetric structure function MRST and MRS(R2), respectively. Figure ^ clearly 
shows that W asymmetry measurements at RHIC could provide an independent determination of d/u. 



6.3 Strange Sea in the Nucleon 

As discussed earher, an interesting consequence of the meson-cloud model is that the s and s distri- 
butions in the proton could have very different shapes, even though the net amount of strangeness 
in the proton vanishes. By comparing the v and v induced charm production, the CCFR collabora- 
tion found no difference between the s and s distributions ||158|| . More precise future measurements 
would be very helpful. Dimuon production experiments using beams might provide an independent 
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Figure 21: Projected statistical accuracy for a {p + d) /2a{p + p) in a two-month run at RHIC 
using the PHENIX detector ||183|| . The E866 data are also shown. 



determination of the s/s ratio of the proton, provided that our current knowledge on valence-quark 
distributions in kaons is improved. As discussed in Section 5.1, ongoing measurements of Ffg via parity- 
violating electron-nucleon scattering should shed much light on the possible difference between s and s 
distributions. 



6.4 Sea Quark Polarization 

Polarized DY and production in polarized p + p collision are planned at RHIC ||1 74|| and they have 
great potential for providing qualitatively new information about antiquark polarization. At large xp 
region {xp > 0.2), the longitudinal spin asymmetry in the p + p DY process is given by ||175| , |176 | 

A.U 

AH{xi, X2) ^ gi{xi)/Fi{xi) X — (xa), (47) 

u 

where gi{x) is the proton polarized structure function measured in DIS, and A'u(x) is the polarized u 
distribution function. 

Equation ^ shows that u polarization can be determined using polarized DY at RHIC. Additional 
information on the sea-quark polarization can be obtained via production ||177|| . The parity- 
violating nature of W production implies that only one of the two beams need to be polarized. At 



positive xp (along the direction of the polarized beam), one finds ||1 77 



AY^^—{x,), and <-^^(x2), (48) 
u a 

where AY is the single-spin asymmetry for W production. Equation ^ shows that the flavor dependence 
of the sea-quark polarization can be revealed via production at RHIC. A remarkable prediction of 
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Figure 22: Predictions of a{p + p — > W'^X)/a{p + p W~X) as a function of Xi? at -y/i = 
500 GeV. The dashed curve corresponds to the d/ u symmetric MRS SO' structure functions, 
while the solid and dotted curves are for the d/u asymmetric structure function MRST and 
MRS(R2), respectively. 



the chiral quark-soliton model is that the flavor asymmetry of polarized sea-quark, Au{x) — Ad{x), is 
very large | p.78|| . This is in striking contrast to the meson cloud model which predicts very small values 

Future DY and 



for Au{x) — Ad(x) 
these models |181| 



179, 180 



production experiments at RHIC could clearly test 



6.5 Hard Processes Tagged by Leading Particles 

The role of the meson-cloud model in explaining the d/u asymmetry suggests a novel technique to study 
meson substructures without using a meson beam. The idea is that the meson cloud in the nucleon 
could be considered as a virtual target to be probed by various hard processes. Recently at the HERA 
e-p collider, meson structure functions were measured in a hard diffractive process, where forward-going 



neutrons or protons were tagged in coincidence with the DIS events [ [L82|| . Analogous measurements 



could be done at RHIC for p-p collisions [|183|| . In particular, a DY pair in coincidence with a forward- 



going neutron or proton could provide information on the antiquark distributions in pions at small x. 
The underlying picture of the tagged DY process is as follows. A proton fluctuates into a baryon plus a 
meson {p—>-n + 7T^,p—>-A + , for example) and the other proton beam interacts with the meson, 
producing a lepton pair. The remaining baryon moves roughly along the initial proton beam direction 
and could be detected at forward angles. A simulation for the p + p ^ n + + X for the PHENIX 
detector, where the muon pairs are detected in the muon arms and the neutrons are detected by a 
small-angle calorimeter showed that such measurement is quite feasible [|183| . 



Drell-Yan experiment tagged by forward-going baryons at RHIC would provide a direct test of the 
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meson-cloud model. One could also conceive tagging forward-going A or A. The A-tagging is of special 
interest since it can shed light on the strange-quark contents of the proton. Another extension is the 
measurement of double- helicity asymmetry, All, in ncutron-taggcd p -|- p Drell-Yan process. If the 
dominant underlying process is p -|- tt interaction, All is expected to be zero. 

7 Conclusion 

The flavor asymmetry of the nucleon sea has been clearly established by recent DIS and DY experiments. 
The surprisingly large asymmetry between u and d is unexplained by perturbative QCD. We draw the 
following conclusions: 

• The X dependence of dp/up indicates that a d,u symmetric sea dominates at small {x < 0.05) and 
large x {x > 0.3). But for OA < x < 0.2 a large and significant fiavor non-symmetric contribution 
determines the sea distributions. 

• The value of the integrals 

J\dp{x, Q^) - Up{x, Q^)]dx = J\un{x, Q^) - dn{x, Q^)]dx = 0.118 ± 0.012 (49) 

are approximately independent of and are characteristic of the nucleon. 

• The value of the integral and the distribution in x can be reproduced using virtual pions. The 
X distribution favors the virtual pion being emitted by the nucleon as a whole, rather than by a 
constituent quark. 

• Virtual pion emission reduces the spin carried by quarks but is not sufficient to account for the 
observed reduction. 

• The mean square charge in the neutron sea is less than that in the proton sea, because > dn 
and dp > Up. 

The up-down flavor asymmetry in the nucleon sea is as fundamental as the axial vector coupling 
constant, qa, in characterizing the properties of the nucleon. This asymmetry is of a completely non- 
perturbative origin but, differently from qa, could only be observed by measurements at the partonic 
level. Thus it stands as a key property of the nucleon that awaited the discovery of perturbative QCD 
processes and the associated high-energy facilities to uncover it. While the nucleon "spin crisis" involves 
both perturbative (spin carried by gluons) and non-perturbative processes (spin flips due to the emission 
of virtual GBs), the origin of the up-down asymmetry in the nucleon sea is entirely non-perturbative. 
Fortunately, it is most likely that the origin of this asymmetry lies with the presence of virtual isovector 
mesons, mostly pions, in the nucleon sea. The unfortunate aspect of this problem is that 70 years of 
research in nuclear physics and the structure of the nucleon is still unable to provide a quantitative 
answer to this basic question. Probably a clear picture would eventually emerge after some future 
experimental investigation which includes: 

• Measurement of dp/up for x > 0.25 and x < 0.01 to more firmly establish the origin of the fiavor 
asymmetry. 

• Ad and Am in the interval 0.1 < a; < 0.2 needs to be measured. Each must be ~ if the pionic 
explanation of the d, u asymmetry is correct. The chiral soliton model predicts that Ad — Au ~ 
—0.5 in the proton for this x interval. 
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• Semi-inclusive including polarized DIS and Drell-Yan experiments to demonstrate the role of 
mesonic degrees of freedom in the nucleon's makeup. 

• More precise measurements on the s versus s distributions in the nucleon. 

• Experimental identification of instantons. 
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